clear all;
close all;
leo1=imread('pepsi1.jpg');%读入图片
leo2=imread('pepsi2.jpg');
% k=融合系数 w=小波系  
k1=0.5;k2=0.5;w='db4'; tic;
outdoor1=rgb2gray(leo1);
outdoor2=rgb2gray(leo2);
%三层小波分解
[ca11,chd11,cvd11,cdd11]=dwt2(outdoor1,w);
[ca12,chd12,cvd12,cdd12]=dwt2(ca11,w);
[ca13,chd13,cvd13,cdd13]=dwt2(ca12,w);
[ca21,chd21,cvd21,cdd21]=dwt2(outdoor2,w);
[ca22,chd22,cvd22,cdd22]=dwt2(ca21,w);
[ca23,chd23,cvd23,cdd23]=dwt2(ca22,w);
%求边缘图像
%矩阵融合
chd3=maxabs(chd13,chd23);
cvd3=maxabs(cvd13,cvd23);
cdd3=maxabs(cdd13,cdd23);
chd2=maxabs(chd12,chd22);
cvd2=maxabs(cvd12,cvd22);
cdd2=maxabs(cdd12,cdd22);
chd1=maxabs(chd11,chd21);
cvd1=maxabs(cvd11,cvd21);
cdd1=maxabs(cdd11,cdd21);
ca3=k1*ca13+k2*ca23;
%逆小波变换
L2=size(chd2);L1=size(chd1);
ca2=idwt2(ca3,chd3,cvd3,cdd3,w);
ca1=idwt2(ca2(1:L2(1),1:L2(2)),chd2,cvd2,cdd2,w);
I=idwt2(ca1(1:L1(1),1:L1(2)),chd1,cvd1,cdd1,w);
toc;

figure(1);
imshow(leo1);
figure(2);
imshow(leo2);
figure(3);
imshow(uint8(I));
title('系数绝对值取最大');